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METHOD AND SYSTEM FOR QUANTITATIVE IMAGE CORRECTION FOR 
OPTICAL COHERENCE TOMOGRAPHY 



CROSS-REFERENCE TO RELATED APPLICATIONS 

This application claims priority under 35 U.S.C. §119 from Provisional 
Application Sen No. 60/395,597 filed July 12, 2002, the entire disclosure of which is 
incorporated herein by reference. 

TECHNICAL FIELD 

The present invention relates generally to the field of optical coherence 
tomography and, more particularly, to a method and system for quantitative image 
correction of optical coherence tomography images of layered media. 



BACKGROUND 

Optical coherence tomography (OCT) is a technology that allows for non-invasive 
micron-scale resolution imaging in living biological tissues. Recent OCT research has 
focused on developing instrumentation appropriate for imaging in clinical settings (e.g., in 

20 ophthalmology, dermatology and gastroenterology), on resolution improvements, real- 
time imaging, and on functional imaging, such as in color Doppler OCT. 

Current-generation real-time OCT systems typically employ depth-priority 
scanning, with the axial scan implemented using a rapid-scan optical delay (RSOD) in the 
reference arm. The rapid axial scan is readily implemented using resonant scanners. 

25 However, the resulting sinusoidally varying delay axially distorts the resulting OCT 
imagines. In addition, the use of non-telecentric scan patterns is often necessitated by 
non-planar sample configurations (e.g., imagining the convex surface of the cornea or the 
concave surface of a hollow organ or tract). 

One major impediment to the use of OCT for quantitative morphological imaging 

30 is image distortions that may occur due to several mechanisms, including nonlinearities in 
the reference or sample scan mechanisms, non-telecentric (diverging or converging) scan 
geometries, and the refraction of probe light in the sample. Non-contact imaging, one of 
the primary advantages of OCT, also leads to significant image distortions due to 
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refraction of the probe beam at the interface between air and smooth surfaces, such as the 
cornea, or liquid accumulations in internal organs. Image distortions due to refraction 
may also occur at internal smooth tissue index boundaries, such as the cornea-aqueous 
interface in the eye. 

5 Accordingly, a need exists for an improved method and system for quantitative 

imagine correction of OCT images, which overcome the above-referenced problems and 
others. 

SUMMARY OF THE INVENTION 

10 According to one aspect of the invention, the invention is directed to a method of 

correcting optical coherence tomography (OCT) data obtained from a layered media 
having at least one interface. The method can include identifying the at least one interface 
from the obtained OCT data and correcting the OCT data for distortion at the at least one 
interface. 

1 5 According to another aspect of the invention, the invention is directed to a 

quantitative image correction method for optical coherence tomography (OCT). The 
method can include correcting for external distortions caused by scan geometry and 
correcting for intrinsic distortions within a sample. 

According to another aspect of the invention, the invention is directed to a non- 
20 invasive system for imaging an anterior portion of an eye. The system can include an 

optical coherence tomography (OCT) data acquisition system and an OCT data correction 
processor, which (i) receives OCT data from the OCT data acquisition system, (ii) 
automatically segments anatomical structures in the anterior portion of the eye to detect at 
least one interface, and (iii) corrects for refraction effects at the at least one detected 
25 interface. 

According to another aspect of the invention, the invention is directed to a non- 
invasive method for imaging an anterior portion of an eye. The method can include 
obtaining optical coherence tomography (OCT) data from the eye. From the obtained 
OCT data, a position of (i) the epithelium, (ii) the endothelium and (iii) the iris can be 
30 determined. Image data distortions caused by at least one of (i) a first interface including 
the epithelium and (ii) a second interface including the endothelium can be corrected. 
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BRIEF DESCRIPTION OF DRAWINGS 

These and further features of the present invention will be apparent with reference 
to the following description and drawings, wherein: 
5 FIG. 1 is a schematic illustration of a non-invasive imaging system for use in 

accordance with the present invention; 

FIG. 2 is a schematic illustration of an exemplary optical coherence tomography 
(OCT) imaging system for use in accordance with the present invention; 

FIG. 3 is a schematic illustration of an exemplary lateral scanning hand held probe 
1 0 for use in accordance with the present invention; 

FIG. 4 is a flow chart illustrating a method for correcting extrinsic and intrinsic 
distortions in OCT data in accordance with the present invention; 

FIG. 5 is a schematic illustration of a mapping array for use in accordance with the 
present invention; 

15 FIG. 6A is a plot schematically illustrating correction of non-linear axial scanning; 

FIG. 6B is a schematic illustration of a coordinate system for raw OCT image 

data; 

FIG. 6C is an exemplary schematic target coordinate system in a homogeneous 
medium; 

20 FIG. 6D is an exemplary schematic illustration of refractive interfaces in a sample 

altering ray paths; 

FIG. 7 is a flow chart illustrating a method of identifying interfaces within the 
sample in accordance with the present invention; 

FIG. 8 illustrates an exemplary series of OCT images corrected in accordance with 
25 the present invention; 

FIG. 9 illustrates an exemplary image to be used for diagnosing risk of acute 
angle-closure glaucoma 

FIG. 10 illustrates an exemplary image to be used for exact fitting of an 
intraocular lens; 
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FIG. 1 1 illustrates exemplary images of the anterior segment of the eye, including 
the epithelium, endothelium, iris, lens and angular portion, corrected in accordance with 
present invention; and 

FIG. 12 illustrates exemplary partial images of the anterior segment corrected in 
5 accordance with the present invention. 

DISCLOSURE OF INVENTION 

In the detailed description that follows, corresponding components have been 
given the same reference numerals regardless of whether they are shown in different 

1 0 embodiments of the present invention. To illustrate the present invention in a clear and 
concise manner, the drawings may not necessarily be to scale and certain features may be 
shown in somewhat schematic form. 

With reference to FIG. 1, a non-invasive optical imaging system is provided. In 
one embodiment, the system includes an optical coherence tomography (OCT) system 10. 

1 5 Raw or otherwise uncorrected OCT data is transmitted to a data correction processor 12. 
As is discussed more fully below, the data correction processor 12 can correct for both 
external distortions (e.g., distortions relating to the scan geometry) and intrinsic 
distortions within the sample (e.g., refraction). The corrected or target OCT data can then 
be presented on a human-readable display 14 for use in determining and/or evaluating 

20 diagnostic parameters. 

FIG. 2 illustrates an exemplary OCT data acquisition system 10, which can be 
used in connection with the present invention. In one embodiment, the OCT data 
acquisition system 10 includes a source 20, such as a super-luminescent diode (SLD) 
source. The optical radiation or light source 20 can illuminate an OCT interferometer 22, 

25 which can include a beam splitter 24 (such as a fused-taper 50/50 fiber coupler). The 
beam splitter 24 separates the optical radiation received from the source 20 into two 
beams. It is to be appreciated that the beam splitter could be other then a 50/50 or 
balanced fiber coupler, such as an unbalanced fiber coupler (e.g., a/(l-cc)). Optical 
radiation can be transmitted to a reference arm 28, including a reference element 26, such 

30 as a translating or rotating reference mirror. The other beam can be transmitted to a 
sample arm 30, optionally including a beam-steering mirror 30, to focus the combined 
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optical radiation on a sample 32. The reflected light received by the beam splitter 24, 
back from both the reference arm 28 and sample arm 30 is combined and transmitted to 
one or more detectors 40, such as photoreceivers or photodetectors. 

It is to be appreciated that a variety of OCT data acquisitions and data acquisition 
5 devices can be employed without departing from the scope of the present invention. For 
example, FIG. 3, schematically illustrates an exemplary lateral scanning handheld probe 
with a diversion scan. Such an OCT system can include a Fourier-domain rapid-scanning 
optical delay (RSOD) in the reference arm incorporating a resonant scanner oscillating at 
2kHz. For the sake of illustration it is to be appreciated that the following abbreviations 
10 are used in FIG. 3; CL - Collimation Lens, SM - Scanning Mirror, RL - Relay Lens, OL - 
Objective Lens, CG - Calibration Grid, f xx - Focal Length of Lens XX. 

With reference now to FIG. 4, a method for correcting distortions in optical 
coherence tomography (OCT) data is provided. The method begins at step 50 with the 
correction of uncorrected or raw OCT data. It is to be appreciated that, in accordance 
15 with one embodiment, the method finds applicability to correcting both external 

distortions (such as those caused by scan geometry) and intrinsic distortions within the 
sample (such as those found in layered media or samples). At step 52, the raw OCT data 
can be corrected for non-linear scanning distortions. Such distortions can include non- 
linear axial scanning distortions, non-linear lateral scanning distortions and/or non- 
20 telecentric scanning distortions. 

In one embodiment, image transformations can be implemented using either 
forward or backward mapping approaches. In the raw (uncorrected) OCT image, x' and 
y' can be defined to denote the coordinates across and along A-scans (single depth scans), 
x and y are the corresponding coordinates in the target (corrected) image (see FIG. 5). 
25 Pixel positions are abbreviated as P' = (x',y') and P = (x,y). A forward mapping approach 
(P = f(P')) calculates the target position for a given raw image data point. Since the target 
position P will most likely fall between target pixels, algorithms can be applied to 
distribute the value of the raw image pixel onto several neighboring pixels. Since the 
density of this distribution is in general nonhomogeneous, the brightness of the target 
30 image must be renormalized, which can be a significant computational expense. 
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A backward mapping approach (P' = F(P)) avoids this disadvantage by mapping 
each target pixel to a location in the raw image, then uses simple interpolations between 
surrounding pixels in the raw image to obtain the target pixel value. Furthermore, the 
number of transformations is at the optimum of one per target pixel, another advantage 
compared with forward mapping. If the same backward transformation is applied to all 
images, it can be implemented with lookup table (also referred to as a pointer array or 
mapping array) to achieve real-time imaging. 

Transformations can be cascaded, using the output of one transformation as the 
input for the next. 

With reference to FIG. 6 and continued reference to FIG. 4, in one embodiment, 
the step of correcting non-linear scanning distortions 52 can include correcting distortions 
due to non-linear axial scanning. The motion of the resonant scanner can be well 
described as sinusoidal, simplifying the backward transformation. In a standard high- 
speed OCT system the carrier is generated by the scan mirror motion. Therefore, the 
nonlinear motion also limits the usable duty cycle r| if no tracking filters are employed. 
Within this duty cycle, the scan depth d is acquired. With this knowledge, the forward 
transformation for the axial scan coordinate y can be written as (FIG. 6A) 




The loial scan depth A can be calculated from 





O) 



Therefore, the backward transformation can he written as 



(4) 
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since the horizontal scan is linear. The transformations to correct non-telecentric 
scanning 

and refraction can use the target image of this transformation as the raw image or the 
transformations can be analytically cascaded. 
5 To correct for geometric image distortions due to non-telecentric scanning, the 

coordinate systems for the raw and target images are defined in terms of the scan 
geometry (FIG 6B, 6C). It can be assumed that the OCT image is formed using a sample 
arm probe that incorporates a rotating mirror scanner [2]. The origin for both systems is in 
the center of the image, x' is linear with acquisition time and thereby with the lateral 

10 scanning mirror movement. The center of the field of view (FOV, having width w and 
depth d) is a distance D from the scanning pivot (SP) or an image of the sample pivot 
(ISP), depending upon the internal details of the sample arm probe. For positive or 
negative D, scans diverge or converge, respectively. For telecentric scans, D approaches 
infinity, but for any real system, this can always be approximated by a large D. 

1 5 The target pixel position P can also be defined in polar coordinates (<t>, L), with 

the scanning angle <t> in the FOV and the distance L from ISP to P (FIG. 6B, 6C). For a 
homogeneous sample without refraction (denoted by subscript h), <t> and L are given by 



The scanning angle to reach the extreme of the FOV at the center plane is given by 

20 <t>max = <t>h. In the rectangular array of acquired data, the scan angle <t> f in the FOV is 
linear with position x': <t>'(x\ y') = 2x'0' m ax/w=X , D, while the distance between ISP 
and P' is L'(x\ y')=:D-y\ Since 4>= O', L = L\ and O max = *max, the complete 
backward transformations are given by: 



<fh ( x, y J = a rctan U f(D - y J j * 



and 





25 
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It is to be appreciated that the step of correcting non-linear scanning distortions 52 
may be performed at various times during the correction process or not at all. 

Referring again to FIG. 4, at step 54, the raw OCT data or the OCT data corrected 
for non-linear scanning distortions can be processed to identify any interfaces present 
5 within the sample. As shown in FIG. 7, identifying the interfaces 54 can include 

thresholding the image data 60. In one embodiment, the image data can be thresholded 
four standard deviations above the noise floor. Alternatively, other appropriate 
thresholding can be performed in order to remove most "salt and pepper noise" from the 
image data. Following the thresholding step 60, the raw image data can be transformed 

10 into a binary image 62 using known techniques. 

Once the image data has been transformed into binary image data, the image data 
can be searched sequentially 64, e.g., on a column-by-column basis. In searching the 
image data 64, positive and negative transitions are sought as indicators of the existence 
of upper and lower interface boundaries. For example, a column of data can be searched 

1 5 wherein each pixel is compared to a threshold indicative of the noise floor for the image 
data. Once a pixel having a value higher than the noise floor is found, that pixel is 
assigned as a potential interface boundary. As the columns of data are searched, interface 
data points are assigned 66 based upon predetermined anatomical rules. For example, in 
the case of imaging the anterior segment of the eye, certain anatomical rules are known. 

20 These rules include: (a) the epithelial and endothelial interfaces are smooth curves that 
are, to a first approximation, parabolic; (b) the epithelial and endothelial interfaces are 
both concave inward; (c) the endothelial interface lies at least 0.3mm below the epithelial 
interface; (d) the iris is disposed below the endothelial interface; (e) an anterior chamber 
is greater than about 0.1mm deep, said endothelial interface and said iris meeting in the 

25 anterior chamber angles; and (f) a lens capsule is apparent within the iris lumen. Of 
course, other anatomical rules could be applied depending on the particular imaging 
application. 

Therefore, in the above ocular example, the columns of data can be searched and 
the first data point determined to be above the noise floor is assigned as a potential point 
30 indicative of the epithelium or epithelial interface (i.e., the upper surface of the cornea). 
Once the potential interfaces are located as set forth in step 66, a first predetermined 
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geometric fit is applied 68 to the assigned interface points of the interfaces. In the ocular 
example, a parabolic fit is applied to the interface points assigned to the epithelium. Of 
course, other geometric fits could be applied depending upon the predetermined 
anatomical rules associated with the particular imaging application. 
5 After the first geometric fitting step 68, a plurality of rejection rules are applied 70 

to the assigned interface points. These rejection rules can be based on predetermined 
anatomical rules and knowledge in connection with the particular imaging application. 
For example, in the case of imaging the anterior segment of the eye, certain rejection rules 
can be applied. These rules include: (a) reject interface data points more than a 

1 0 predetermined distance apart from the first parabolic fit; (b) reject data points without 
immediate left and right neighboring data points; (c) reject data points that, along with 
neighboring data points, are within a top 0.1mm of the overall image; (d) reject iris points 
about estimated endothelial and epithelial interfaces and below estimated iris; and (e) for 
iris data points, replace vertical positions by the median of neighborhood data points. Of 

15 course, other rejection rules could be applied depending upon the predetermined 

anatomical rules associated with the particular imaging application. It is to be appreciated 
that the above rejection rules could be applied through one or more iterations in which 
tolerances for rejected points are steadily and progressively decreased. 

After the plurality of rejection rules are applied 70, a second predetermined 

20 geometric fit 72 can be applied to the remaining assigned interface points to provide a 
second estimate of the interface or corresponding interface structure. In the ocular 
example discussed herein, the second fit could be another parabolic fit or, alternatively, a 
quadrabolic fit. This second fit and/or estimate 72 should identify the present interface 
74. 

25 Returning to the example of imaging the anterior segment of the eye, the above- 

described methodology for identifying interfaces (illustrated in FIG. 7) would also be 
repeated for the endothelium or endothelial interface, which comprises the lower 
boundary of the cornea. 

By way of example, FIG. 8 provides exemplary images illustrating the 

30 methodology described in FIG. 4 and FIG. 7 in connection with an ocular example 

described above. For example, raw OCT data (upper left image) is corrected for non- 
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linear scanning distortions (upper right image). The image is thresholded (middle left 
image) and binary morphed (middle right image). The initial points are identified and 
assigned (bottom left image) and then certain points are rejected (bottom right image). 

In addition, it is to be appreciated that the interface identifying (also referred to as 
5 segmentation) methodology discussed above and illustrated in FIG. 7, can provide useful 
information about the anterior segment and chamber of the eye. This information can be 
employed in a number of diagnostic applications, including, but not limited to 
determination of angle parameters related to the intersection of the endothelium and the 
iris, which can be useful in evaluating the risk of acute angle-closure glaucoma, and the 

10 exact fitting of intraocular lenses. Appositional angle closure is a risk factor for 

progressive trabecular damage, elevated intraocular pressure and acute angle-closure 
glaucoma. Phakic intraocular lenses (IOL) are a relatively new class of devices implanted 
to correct severe myopia and hyperopia. The exact dimensions of the anterior chamber 
(AC) must be known to determine the lens size for a correct fit. In accordance with the 

15 methodologies discussed herein, OCT is capable of rapidly imaging the complete AC with 
a resolution of about 10 microns. 

FIG. 9 and FIG. 10 provide exemplary images relating to diagnosing risk of acute angle- 
closure glaucoma (FIG. 9) and exact fitting of intraocular lenses (FIG. 10). 

Referring again to FIG. 4, once the interfaces within the sample (e.g., the epithelial 

20 interface and the endothelial interface) have been identified 54, this information can be 
used to provide correction for the internal distortion 56 (e.g., refraction) caused by these 
interfaces. While the present invention is most concerned with correction of refraction 
effects caused by interfaces within the sample, it is to be appreciated that the present 
invention may be amenable to correction of other intrinsic distortions, such as diffraction 

25 and the like. 

In order to correct for refraction, the refractive indices of sample layers must be 
known and the interfaces between them must be identified (FIG. 6D). Then, fundamental 
principles can be applied to correct for beam deflection and to transform optical into 
physical pathlengths. In the following, it is assumed k layers of different media in the 
30 image area, having indices of refraction n t to n k , with the probe beam entering the sample 
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from the top of the image in air (no = 1). The interfaces I k (x) between layers of differing 
index are assumed to be smooth functions. 

In one embodiment, a forward transformation for refraction correction could use 
Snell's law to calculate the target pixel for each raw data pixel, by propagating the 
5 incident beam sequentially through the sample layers and keeping track of optical path 
length until the correct distance is reached. However, the raw image is distorted by the 
scan geometry and refraction on previous interfaces. Therefore, defining the normal on 
the interface and the incident angle to apply Snell's law becomes very complicated in raw 
image coordinates. For the backward transformation, the interfaces can be defined 

10 distortion-free in target coordinates, avoiding this difficulty. However, for the latter 
transformation, Snell's law cannot be applied since the route of the beam through the 
sample is not known a priori. A solution can be obtained, however, by applying Fermat's 
principle, which states that light rays will propagate such that the optical path between 
source and target locations is minimized. Assuming the sample beam passes though the 

15 points Pi at the interfaces of the layers to reach P (FIG. 6D), the total optical pathlength is 
the sum of the distance from ISP to Pi (Eq. (7)) and the optical path between subsequent 
Pi's to P: 



By varying the location of Pj along the interfaces, a minimum of L can be found, 
20 satisfying 

Fermat's principle. Assuming the paths of the different probe beams do not cross within 
the 

FOV, a unique solution exists for the Pj. After iteratively solving for the values of Pi, the 
complete back transformation can be written as 




(10) 



x-=F K (P 11 ... f P tf P) = F lfc (P 1 )- 



y , = F y (P 1 ....,P Jk .P) = D-UP 1 ....,P*.P)- 
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Referring again to FIG. 4, once the OCT data has been corrected for refraction 56 
or other distortions at the interfaces, the corrected image can be displayed. FIG. 1 1 
provides an exemplary image of the anterior chamber of the eye, including the cornea 
(defined by the epithelium and endothelium), the iris, the angular chambers and the lens. 
As described above, the position of the detected epithelium and endothelium of the cornea 
were used for refraction correction. The angle portions of the anterior chamber can be 
detected and refined using a dynamic programming technique. Because the mechanics of 
dynamic programming are known to skilled artisans, the details will not be discussed 
herein. 

As shown in FIG. 12, the present invention is insensitive to partial views such that 
only part of the anterior segment need be imaged (e.g., with a handheld scanner). These 
partial views are processed in the same manner described more fully above. 

Although, particular embodiments of the invention have been described in detail, 
it is understood that the invention is not limited correspondingly in scope, but includes all 
changes, modifications, and equivalents coming within the spirit and terms of the claims 
appended hereto. In addition, it is to be appreciated that features shown and described 
with respect to a given embodiment may also be used in conjunction with other 
embodiments. 
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